A simple model of price formation 
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A simple Ising spin model which can describe the mechanism of price formation in financial markets 
is proposed. In contrast to other agent-based models, the influence does not flow inward from the 
surrounding neighbors to the center site, but spreads outward from the center to the neighbors. The 
model thus describes the spread of opinions among traders. It is shown via standard Monte Carlo 
simulations that very simple rules lead to dynamics that duplicate those of asset prices. 
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The Ising spin system is one of the most frequently 
used models of statistical mechanics. Its simplicity (bi- 
nary variables) makes it appealing to researchers from 
other branches of science including biology B , sociology 
m and economy |^-||. It is rather obvious that Ising- 
type models cannot explain origins of very complicated 
phenomena observed in complex systems. However, it 
is believed that these kind of models can describe some 
universal behavior. 

Recently, an Ising spin model which can describe the 
mechanism of making a decision in a closed community 
was proposed Q . In spite of simple rules the model ex- 
hibited complicated dynamics in one and two [ 10 1 di- 
mensions. In contrast to usual majority rules [11], in 
this model the influence was spreading outward from the 
center. This idea seemed appealing and we adapted it 
to model financial markets. We introduced new dynamic 
rules describing the behavior of two types of market play- 
ers: trend followers and fundamentalists. The obtained 
results were astonishing - the properties of simulated 
price trajectories duplicated those of analyzed historic 
data sets. Three simple rules led to a fat-tailed distri- 
bution of returns, long-term dependence in volatility and 
no dependence in returns themselves. 

We strongly believe that this simple and parameter free 
model is a good first approximation of a number of real 
financial markets. But before we introduce our model we 
review some of the stylized facts about price formation 
in the financial markets. 

Stylized facts. - Adequate analysis of financial 
data relies on an explicit definition of the variables un- 
der study. Among others these include the price and the 
change of price. In our studies the price Xt is the daily 
closing price for a given asset. The change of price Tt 
at time t is defined as r t — logxt+i — logxt- In fact, 
this is the change of the logarithmic price and is often 
referred to as return. The change of price, rather than 
the price itself, is the variable of interest for traders (and 
researchers as well). 

Fat-tailed distribution of returns. - The variety of 
opinions about the distributions of asset returns and their 
generating processes is wide. Some authors claim the 
distributions to be close to Paretian stable |12] , some to 



generalized hyperbolic [^3|,^4| , and some reject any single 
distribution 

Instead of looking at the central part of the distribu- 
tion, an alternative way is to look at the tails. Most 
types of distributions can be classified into three cate- 
gories: 1° thin-tailed - for which all moments exist and 
whose density function decays exponentially in the tails; 
2° fat-tailed - whose density function decays in a power- 
law fashion; 3° bounded - which have no tails. 

Virtually all quantitative analysts suggest that asset 
returns fall into the second category. If we plot returns 
against time we can notice many more outlying (away 
from the mean) observations than for white noise. This 
phenomenon can be seen even better on normal proba- 
bility plots, where the cumulative distribution function 
(CDF) is drawn on the scale of the cumulative Gaus- 
sian distribution function. Normal distributions have the 
form of a straight line in this representation, which is ap- 
proximately the case for the distribution of weekly or 
monthly returns. However, distributions of daily and 
higher- frequency returns are distinctly fat-tailed jl7||. 
This can be easily seen in the top panels of Fig. 1, where 
daily returns of the DJIA index for the period Jan. 2nd, 
1990 - Dec. 30th, 1999, are presented. 

Clustering and dependence. — Despite the wishes of 
many researchers asset returns cannot be modeled ad- 
equately by series of iid (independent and identically 
distributed) realizations of a random variable described 
by a certain fat-tailed distribution. This is caused by 
the fact that financial time series depend on the evo- 
lution of a large number of strongly interacting sys- 
tems and belong to the class of complex evolving sys- 
tems. As a result, if we plot returns against time 
we can observe the non-stationarity (heteroscedastic- 
ity) of the process in the form of clusters, i.e. peri- 
ods during which the volatility (measured by standard 
deviation or the equivalent I 1 norm [Q) of the pro- 
cess is much higher than usual, see the top-left panel 
of Fig. 1. Thus it is natural to expect dependence 
in asset returns. Fortunately, there are many meth- 
ods to quantify dependence. The direct method con- 
sists in plotting the autocorrelation function: acf (r, k) = 

Etlfe+i( r * ~ f )( r t-fc - r)/J2^ =1 (n - r) 2 , where N is the 
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sample length and f = j? 2tLi r ti f° r different time lags 
k. For most financial data autocorrelation of returns dies 
out (or more precisely: falls into the confidence interval 
of Gaussian random walk) after at most a few days and 
long-term autocorrelations are found only for squared or 
absolute value of returns [is] |2^| , see bottom panels of 
Fig. 1. Recall that for Brownian motion - the classical 
model of price fluctuations jy^^l) - autocorrelations of 
Tt-, r\ and \rt\ are not significant for lags greater or equal 
to one. 
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FIG. 1. Daily returns of the Dow Jones Industrial Average 
index during the last decade (top left), normal probability plot 
of DJIA returns (top right), lagged autocorrelation function 
of DJIA daily returns (bottom left), and lagged autocorrela- 
tion function of absolute value of DJIA daily returns (bottom 
right). Dashed horizontal lines represent the 95% confidence 
interval of a Gaussian random walk. 

Another way to examine the dependence structure is 
the power spectrum analysis, also known as the "fre- 
quency domain analysis". One of the most often used 
techniques was proposed by Geweke and Porter-Hudak 
p2| (GPH) and is based on observations of the slope of 
the spectral density function of a fractionally integrated 
series around the angular frequency lu = 0. A simple 
linear regression of the periodogram /„ (a sample ana- 
logue of the spectral density) at low Fourier frequencies 
u k : \og{I n (u> k )} = a-dlog{4sin 2 (w fc /2)} + e fe yields the 
differencing parameter d = H — 0.5 through the rela- 
tion d = d. The GPH estimate of the Hurst exponent 
H has well known asymptotic properties and allows for 
construction of confidence intervals p^ , p^ |. 

Yet another method is the Hurst R/S analysis [ p4|j25| 
or its " youn ger sister" - the Detrended Fluctuation Anal- 
ysis (DFA) |26|. Both methods are based on a similar 
algorithm, which begins with dividing the time series 
(of returns) into boxes of equal length and normalizing 
the data in each box by subtracting the sample mean 
(R/S) or a linear trend (DFA). Next some sort of volatil- 
ity statistics is calculated (rescaled range or mean square 



fluctuation, respectively) and plotted against box size on 
a double-logarithmic paper. Linear regression yields the 
Hurst exponent H , whose value can lie in one of the three 
regimes: 1° H > 0.5 - persistent time series (strict long 
memory), 2° H = 0.5 - random walk or a short-memory 
process, 3° H < 0.5 - anti-persistent (or mean-reverting) 
time series. Unfortunately, no asymptotic distribution 
theory has been derived for the R/S and DFA statistics 
so far. However, it is possible to estimate confidence in- 
tervals based on Monte Carlo simulations p3[ . 

The model. — Recently a simple model for opinion 
evolution in a closed community was proposed p|,|i"c|] . In 
this model (called USDF) the community is represented 
by a horizontal chain of Ising spins which are either up or 
down. A pair of parallel neighbors forces its two neigh- 
bors to have the same orientation (in random sequential 
updating), while for an antiparallel pair, the left neigh- 
bor takes the orientation of the right part of the pair, 
and the right neighbor follows the left part of the pair. 

In contrast to usual majority rules in the USDF 
model the influence does not flow inward from the sur- 
rounding neighbors to the center site, but spreads out- 
ward from the center to the neighbors. The model thus 
describes the spread of opinions. The dynamic rules lead 
to three steady states: two ferromagnetic (all spins up or 
all spins down) and one antiferromagnetic (an up-spin is 
followed by a down-spin, which is again followed by an 
up-spin, etc.). 

In this paper we modify the model to simulate price 
formation in a financial market. The spins are interpreted 
as market participants' attitude. An up-spin (Si = 1) 
represents a trader who is bullish and places buy orders, 
whereas a down-spin (Si = — 1) represents a trader who 
is bearish and places sell orders. In our model the first 
dynamic rule of the USDF model remains unchanged, 
i.e. if SiSi+i = 1 then Si-i and Si+2 take the direction 
of the pair (i,i+l). This can be justified by the fact 
that a lot of market participants are trend followers and 
place their orders on the basis of a local guru's opinion. 
However, the second dynamic rule of the USDF model 
has to be changed to incorporate the fact that the absence 
of a local guru (two neighboring spins are in different 
directions) causes market participants to act randomly 
rather than make the opposite decision to his neighbor: 
if SiSi+i = — 1 then Si-± and Si+2 take one of the two 
directions at random. 

Such a model has two stable states (both ferromag- 
netic), which is not very realistic for a financial market. 
Fortunately, trend followers are not the only participants 
of the market 27 1. There are also fundamentalists - play- 
ers that know much more about the system and have a 
strategy (or perhaps we should call them "rationalists"). 
To make things simple, in our model we introduce one 
fundamentalist. He knows exactly what is the current 
difference between demand and supply in the whole sys- 
tem. If supply is greater than demand he places buy 
orders, if lower - sell orders. 

It is not clear a priori how to define the price in a 
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market. The only obvious requirement is, that the price 
should go up, when there is more demand than supply, 
and vice versa. For simplicity, we define the price Xt in 
our model as the normalized difference between demand 
and supply (magnetization): x t — f Eili^Wi where 
N is the system size. Note that x t G [—1, 1], so \x t \ can be 
treated as probability. Now we can formulate the third 
rule of our model: the fundamentalist will buy (i.e. take 
value 1) at time t with probability \x t \ if x t < and sell 
(i.e. take value —1) with probability \x t \ if x t > 0. 

The third rule means that if the system will be close 
to the stable state " all up" , the fundamentalist will place 
sell orders with probability close to one (in the limiting 
state exactly with probability one) and start to reverse 
the system. So the price will start to fall. On the con- 
trary, when r will be close to —1, the fundamentalist will 
place buy orders (take the value 1) and the price will 
start to grow. This means that ferromagnetic states will 
not be stable states anymore. 

Results. To investigate our model we perform a 
standard Monte Carlo simulation with random updating. 
We consider a chain of N Ising spins with free boundary 
conditions. We were usually taking N — 1000, but we 
have done simulations for N = 10000 as well. We start 
from a totally random initial state, i.e. to each site of 
the chain we assign an arrow with a randomly chosen 
direction: up or down (Ising spin). 
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FIG. 2. A typical path of the simulated price process Xt 
and the USD/DEM exchange rate, respectively. 

In our simulations each Monte Carlo step (MCS) rep- 
resents one trading hour; eight steps constitute one trad- 
ing day. We typically simulate 20000 MCSs, which cor- 
responds to 2500 trading days or roughly 10 years. We 
chose such a period of time, because in this paper we com- 
pare the empirical results with historical data sets (two 
FX rates and two stock indices) of about the same size: 
2245 daily quotations of the dollar-mark (USD/DEM) ex- 
change rate for the period Aug. 9th, 1990 - Aug. 20th, 
1999 (see Fig. 2 and Table 1), 2809 daily quotations 
of the yen-dollar (JPY/USD) exchange rate for the pe- 
riod Jan. 2nd, 1990 - Feb. 28th, 2001 (see Table 1), 
2527 daily quotations of the Dow Jones Industrial Aver- 
age (DJIA) index for the period Jan. 2nd, 1990 - Dec. 
30th, 1999 (see Fig. 1), and 1561 daily quotations of the 
WIG20 Warsaw Stock Exchange index (based on 20 blue 
chip stocks from the Polish capital market) for the period 
Jan. 2nd, 1995 - Mar. 30th, 2001 (see Table 1). 



The returns r t are obtained from the simulated price 
curve x t (see Fig. 2) after it is shifted (incremented by 
one) to make it positive. Alternatively we could have de- 
fined the up-spin to be equal to two and the down-spin 

to zero. However, this would have made the calcula- 
tions more difficult and the description of the model less 
appealing. 

In Figure 3 we present daily returns and normal proba- 
bility plots for the simulated (left panels) and USD /DEM 
exchange rate (right panels) time series. Without prior 
knowledge as to the magnitude of historical returns it is 
impossible to judge which process is real and which is a 
fraud. The same is true for the simulated and the DJIA 
returns of Fig. 1. 
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FIG. 3. Returns of the simulated price process x t and 
daily returns of the USD/DEM exchange rate during the last 
decade, respectively {top panels). Normal probability plots of 
rt and USD/DEM returns, respectively, clearly showing the 
fat tails of price returns distributions {bottom panels). 

In Figure 4 we present the lagged autocorrelation study 
for the simulated (left panels) and USD/DEM exchange 
rate (right panels) time series. Again the properties of 
the simulated price process duplicate those of historical 
data. The lagged autocorrelation of returns falls into 
the confidence interval of Gaussian random walk imme- 
diately, like for the dollar-mark exchange rate. The same 
plot for the DJIA returns shows a bit more dependence. 

This can be seen also in Table 1, where values and 
significance of the R/S, DFA and GPH statistics for the 
simulated and four historical data sets are presented. In 
all but one case (DFA for DJIA) the Hurst exponents 
of daily returns are insignificantly different from those 
of white noise. On the other hand, the Hurst exponents 
of the absolute value of daily returns are persistent in 
all cases, with the results being significant even at the 
two-sided 99% level. Moreover, the estimates of H from 
the simulation are indistinguishable from those of real 
market data. 
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FIG. 4. Lagged autocorrelation functions of rt and 
USD/DEM returns, respectively (top panels). Lagged auto- 
correlation functions of absolute value for the same data sets 
(bottom panels). Dashed horizontal lines represent the 95% 
confidence interval of a Gaussian random walk. 

The presented empirical analysis clearly shows that 
three simple rules of our model lead to a fat-tailed dis- 
tribution of returns, long-term dependence in volatility 
and no dependence in returns themselves as observed for 
market data. Thus we may conclude that this simple 
model is a good first approximation of a number of real 
financial markets. 

We gratefully acknowledge critical comments of an 
anonymous referee, which led to a substantial improve- 
ment of the paper. 
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Data 


R/S-AL 


Method 
DFA 


GPH 






Returns 




Simulation 


0.5270 


0.4666 


0.3653 


USD/DEM 


0.5127 


0.5115 


0.6154 


JPY/USD 


0.5353 


0.5303 


0.5790 


DJIA 


0.4585 


0.4195** 


0.3560 


WIG20 


0.5030 


0.4981 


0.4604 




Absolute value of 


returns 


Simulation 


0.8940*** 


0.9335*** 


0.8931*** 


USD/DEM 


0.7751*** 


0.8406*** 


0.8761*** 


JPY/USD 


0.8576*** 


0.9529*** 


0.9287*** 


DJIA 


0.7838*** 


0.9080*** 


0.8357*** 


WIG20 


0.9103*** 


0.9494*** 


0.8262*** 



*, ** and *** denote significance at the (two-sided) 90%, 95% 
and 99% level, respectively. For the R/S-AL and DFA statis- 
tics inference is based on empirical Monte Carlo results of 
Weron J2^], whereas for the GPH statistics - on asymptotic 
distribution of the estimate of H [E2|. 
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